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Abstract 

We review the impact of massive neutrinos on cosmological observables at the linear order. By means of N-body 
simulations we investigate the signatures left by neutrinos on the fully non-linear regime. We present the effects 
induced by massive neutrinos on the matter power spectrum, the halo mass function and on the halo-matter bias in 
massive neutrino cosmologies. We also investigate the clustering of cosmic neutrinos within galaxy clusters. 
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1. Introduction 

The Big Bang theory predicts the existence of a cos¬ 
mic neutrino background (CNB). Since neutrinos de¬ 
coupled from the primordial plasma when they were 
ultra-relativistic, their momenta will follow a Fermi- 
Dirac distribution with a temperature equal to T v (z ) = 
(4/ll) 1/3 r C MB(z), where T C MBfe) is the temperature of 
the cosmic microwave background (CMB). 

At the linear order, massive neutrinos impact on cos¬ 
mology in two different ways O : modifying the matter- 
radiation equality time and slowing down the growth of 
matter perturbations. It can be shown that these two ef¬ 
fects induce a suppression, on small scales, in the matter 
power spectrum with respect to a cosmology with mass¬ 
less neutrinos. The above effects are commonly used to 
put tight constrains on the neutrino masses from cosmo¬ 
logical observables such as the anisotropies in the CMB 
and the spatial distribution of galaxies at low redshift. 

While the signatures left by massive neutrinos at the 
linear level have allowed us to put tight constrains on 
their masses, a huge amount of information, contained 
in non-linear scales, is not used since their theoreti¬ 
cal counterpart has not been fully investigated. There¬ 
fore, in order to extract the maximum information from 
the cosmological observables it is crucial to study the 
impact of massive neutrinos into the fully non-linear 


regime. There are different ways to carry this task out: 
by means of perturbation theory, semi-analytic models 
or N-body simulations. Among those, the most precise 
one is the use of N-body simulations, although it is also 
the most computationally expensive. 

2. Signatures in the non-linear regime 

We have run N-body simulations using the code 
GADGET-III [1; massive neutrinos are simulated as 
additional particles. Our simulation suite comprises 
three different cosmologies: 2 m v = 0-0 eV, 2 = 0-3 

eV and 2 wiv = 0-6 eV. We note that even thought mod¬ 
els with Yu m v > 0.15 eV are ruled out by the latests 
analysis that combine CMB, BAO and Lyman-cr forest 
data (HID, our purpose is to investigate the impact of 
massive neutrinos into non-linear observables as a func¬ 
tion of the neutrino masses; thus, we are justified to con¬ 
sider such cosmological models. 

In Fig. [I] we show the total matter spatial distribution 
as obtained from our simulations, at z = 0, for two dif¬ 
ferent cosmologies: a model with Yj m v = 0.0 eV (top) 
and a model with £ = 0.6 eV (bottom), the left pan¬ 

els show a zoom into the region marked with a purple 
square in the right panels. The differences among the 
two models can be seen directly from visual inspection. 
In particular, it is interesting to note that the model with 
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Figure 1: Spatial distribution of the matter density field in a universe 
with Yj m v - 0-6 eV (top) and massless neutrinos (bottom) taken from 
a slice of 2 h ~ l Mpc of our N-body simulations at z = 0. The big 
panels show the distribution on large scales whereas the small panels 
display a zoom into the regions marked in purple. 


massless neutrinos is a more ’’evolved” version of the 
model with 2 = 0.6 eV neutrinos. This can be seen 

from the bottom-left panel, where two halos are merg¬ 
ing with the big halo located in the middle-upper part of 
that panel. However, for the model with massless neu¬ 
trinos those two halos are further apart, indicating that 
the evolution in that universe is slower, as predicted by 
linear theory. 

The fully non-linear matter power spectrum repre¬ 
sents the most basic observable to investigate the impact 
of massive neutrinos. Several authors have already stud¬ 
ied in detail the neutrino signatures in the matter power 
spectrum on small scales in the non-linear regime by 
means of N-body simulations 0|6). It is found that 
massive neutrinos suppress power on small scales, with 
respect to the equivalent massless neutrino cosmology. 
This suppression is larger than the linear prediction up 
to k ~ 1 /zMpc -1 , whereas on smaller scales the sup¬ 
pression is lower than the linear one. These features can 
be explained using an extension of the halo model as the 
one presented in 0. 

An important observable in which massive neutrinos 
leave their signature is in the halo mass function (HMF), 
i.e. the abundance of dark matter halos of a given mass. 


Theoretical studies have pointed out that the mass func¬ 
tion becomes universal GO once it is expressed in terms 
of the peak height, v = SJcr(M), where 6 C is the crit¬ 
ical overdensity required for collapse ( S c ~ 1.686) and 
<x(M) is the r.m.s. value of the linear density field when 
smoothed with a top-hat filter of mass M. It has been 
a common practice to estimate the HMF in massive 
neutrino cosmologies by computing cr(M) using the to¬ 
tal matter linear power spectrum. Recently, in a series 
of papers (9j ljj, IT ] we have demonstrated using N- 
body simulations that the HMF in massive neutrinos 
cosmologies can be well reproduced if cr(M) is com¬ 
puted with respect to the cold dark matter field, instead 
of the total matter field. Besides, by doing that, the 
HMF becomes universal. The use of an incorrect HMF 
for massive neutrinos cosmologies may bias the inferred 
values of the cosmological parameters as pointed out in 

0 

Massive neutrinos also imprint their signature on the 
bias between the spatial distribution of dark matter halos 
and the one of the underlying matter. It is well known 
that, for massless neutrinos cosmologies, on large, lin¬ 
ear scales, the halo-matter bias should become scale- 
independent; we have confirmed this by using our mass¬ 
less neutrinos N-body simulations. However, in massive 
neutrinos cosmologies we find the halo-matter bias to be 
scale-dependent, even on very large scales. This can be 
understood if we account for the fact that halos trace 
the spatial distribution of cold dark matter, rather than 
all the matter. Our results point out that if the bias is 
defined as the square root of the ratio between the ha¬ 
los and cold dark matter power spectra, then the bias is 
scale-independent and universal. 

Another interesting feature of the CNB is that some 
of the cosmic neutrinos may get trapped within the deep 
gravitational potential wells present in galaxy clusters. 
To see this more clearly let us write down some num¬ 
bers. The velocity dispersion of galaxy clusters, at 
z = 0, is of the order ~ 1000 km/s. On the other hand, 
the mean thermal velocities of cosmic neutrinos is given 
by V v (z) = 160 (^) km/s, where m v is the mass of 
each neutrino eigenstate in eV, i.e. V v (z) = 800 km/s 
for neutrinos with £ Wv = 0-6 eV at z = 0. Therefore, 
neutrinos in the low-velocity tail of the Fermi-Dirac dis¬ 
tribution are expected to cluster within the most massive 
dark matter halos. This clustering of relic neutrinos can 
be seen in Fig. [2] In the left panel of that figure we show 
the overdensity of CDM around one of the most massive 
halos in our simulation with Yj m v — 0.6 eV neutrinos, 
while the right panel instead displays the massive neu¬ 
trinos overdensity. 
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Figure 2: Cold dark matter (left) and neutrino (right) overdensity 
around a massive dark matter halo at z = 0 for a cosmological model 
with Yj m v = 0.6 eV neutrinos. 

As can be seen from Fig. [2j massive neutrinos do 
cluster around massive dark matter halos. There are few 
interesting features to note. The first one is that cos¬ 
mic neutrinos do not cluster as much as the CDM. This 
can be seen from the scale located below the panels of 
Fig. [2] The reason is that only neutrinos from the low- 
velocity tail of the momentum distribution can cluster. 
The second feature is that the neutrino halo it is much 
more spatially extended that its CDM counterpart. This 
can be understood taking into account the large thermal 
velocities associated with neutrinos, which make them 
to have very large orbits. Moreover, those thermal ve¬ 
locities avoid the clustering of neutrinos in the center 
of the dark matter halo, inducing a core in the neutrino 
density profile. 

Authors in G2 studied the shape and amplitude of 
the neutrino density profile using N-body simulations. 
They found that the density profile of the neutrino halos 
can be well fitted by the expression 6 v {r) = (p v (r ) - 
Pv)/Pv = PcK 1 + ( r/r c ) a ), where p c , r c and a are 
free parameters. Moreover, those authors demonstrated 
that that profile is universal, i.e. valid for all neutrino 
masses, for different host dark matter halo masses...etc. 

3. Summary and Conclusions 

The existence of a cosmic neutrino background is 
predicted by our standard cosmological model. The dy¬ 
namics of those relic neutrinos is very different to the 
one of the cold dark matter, since the thermal velocities 
of the former are very large whereas they are assumed 
to be negligible for the latter. The impact of massive 
neutrinos on cosmological observables is very well un¬ 
derstood at the linear order, while the signatures left by 


neutrinos into the fully non-linear regime need sophisti¬ 
cated, and computationally expensive, tools as N-body 
simulations. 

By running large box size N-body simulations that 
incorporate massive neutrinos as additional particles we 
have investigated the impact of massive neutrinos on: 
the matter power spectrum, the halo mass function, 
the halo-matter bias and the neutrino clustering within 
galaxy clusters. 

We find that the halo mass function in a massive neu¬ 
trinos cosmology can be reproduced from the fitting 
functions of the corresponding massless neutrinos cos¬ 
mology, by computing the relevant quantities (such as 
<x(M)) with respect to the CDM field, instead of the to¬ 
tal matter field. By doing so we find that the halo mass 
function becomes universal. 

The halo-matter bias, defined as the square root of 
the ratio between the halos power spectrum to the to¬ 
tal matter power spectrum, exhibits a scale-dependence 
behavior even on linear scales. This is a unique fea¬ 
ture not present in standard massless neutrinos cosmolo¬ 
gies. Our results point out that if we define the bias as 
the square root of the ratio between the halos and cold 
dark matter power spectra, then the bias becomes scale- 
independent on large scales and universal. 

The momenta of the cosmic neutrino background fol¬ 
low a Fermi-Dirac distribution. Thus, neutrinos from 
the low velocity tail of the distribution are expected 
to cluster within the deep gravitational potential wells 
present in galaxy clusters. In our simulations we find 
halos of neutrinos located in the potential wells of mas¬ 
sive dark matter halos. We find the density profile of 
the neutrino halos to be universal. Possible detection 
of such neutrino halos may be achieved through weak 
lensing lfl3l . 
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